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Abstract. We consider the effect of two different environments on the 
performance of the quantum adiabatic search algorithm, a thermal bath at finite 
temperature, and a structured environment similar to the one encountered in 
systems coupled to the electromagnetic field that exists within a photonic crystal. 
While for all the parameter regimes explored here, the algorithm performance is 
worsened by the contact with a thermal environment, the picture appears to be 
different when considering a structured environment. In this case we show that, by 
tuning the environment parameters to certain regimes, the algorithm performance 
can actually be improved with respect to the closed system case. Additionally, the 
relevance of considering the dissipation rates as complex quantities is discussed 
in both cases. More particularly, we find that the imaginary part of the rates can 
not be neglected with the usual argument that it simply amounts to an energy 
shift, and in fact influences crucially the system dynamics. 
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1. Introduction 

The paradigm of Adiabatic Quantum Computation (AQC), introduced in [1], is known 
to be equivalent to standard, circuit based, quantum computation [2] . From the point 
of view of its physical implementation, it is however specially appealing, since it is 
entirely formulated in terms of a controlled dynamics, and moreover it offers some 
inherent robustness to errors [3]. 

In AQC the solution of a computational problem is encoded in the ground state 
of a certain (problem) Hamiltonian. The system is started in the easily preparable 
ground state of an initial Hamiltonian, which is then slowly turned into the problem 
one. The adiabatic theorem [H [S] El [Zl [H] guarantees that, if the variation of the 
Hamiltonian is sufficiently slow, time evolution will drive the system into the solution 
state. The resource measuring the computational cost of the algorithm is thus the total 
evolution time required to guarantee the adiabaticity conditiorlj. For most cases, it 
depends on the inverse squared of the minimum energy gap during the evolution. 

Since quantum systems are not in general completely isolated from their 
environment, errors and dissipation are ubiquitous, and it is fundamental to decide 
how quantum computers and algorithms can be built that achieve their computational 
tasks in spite of inaccurate operations and certain loss of coherence. But AQC gets 
affected by errors in a fundamentally different way than the standard, gate-based 
model of quantum computation. The interaction with the environment may excite 
the system out of its ground state, causing errors, but if the energy scales in the 
environment are much smaller than the minimum gap, the adiabatic evolution will be 
naturally preserved. On the other hand, the dissipation may also result in a modified 
effective Hamiltonian, and thus affect the performance of the AQC [9j [TOl [11] . Even 
when the coupling with the environment is weak, it may produce noticeable effects in 
the performance of the adiabatic algorithm, specially considering the large time scales 
that this requires. 

Looking at the question from a different perspective, the algorithm may be 
implemented within a system inwhich the interaction with the environment is not 
an undesirable effect, but rather a tool that can be actually controlled and tuned 
at will to improve the performance of AQC. Whether the effect of the environment 
is controllable or not, the system in question has to be considered as a quantum 
open system, and its interaction with the environment can be treated within the weak 
coupling limit, provided that this coupling is small enough in comparison to the system 
and the environment time scales. 

In this work we study the effect of different types of noise on a particular adiabatic 
algorithm, namely the adiabatic version of Grover's search. Grover's problem, or 
that of search in an unstructured database [12], is one of the problems for which 
quantum computation has been explicitly shown to exhibit a remarkable speedup over 
the classical one. The adiabatic version of the quantum search algorithm achieves, 
in a closed system, the optimal quadratic speedup with respect to its best classical 
counterpart [13 . But if the system is subject to some dissipative dynamics, the 
performance of the search may depend on the characteristics of the bath [14] . In 
particular, it was pointed out by Amin and collaborators [10. that the presence of a 
thermal bath can in some cases enhance the performance of the adiabatic quantum 
search. 

J A more strict measure of the cost is given by the product of the norm of the Hamiltonian times 
the total time [2]. 
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In this paper, we study the performance of the adiabatic search algorithm 
implemented in a quantum open system. To this end, we adopt the two-level 
approximation, in which only the interaction of the environment with the two levels 
involved in the quantum search algorithm is considered 114] . 

We will consider two situations. 

First, we analyze the case in which the system is coupled to a thermal 
environment. In particular, we consider an ohmic environment as appearing in |15j . 
which is valid to describe cases in which the spectral density is a smooth function in 
the frequency space. 

Second, we consider an ideal situation in which the adiabatic search algorithm 
is implemented within an atom lattice. In that case, we study the model proposed 
in [TB], where the atoms in the lattice are coupled in a controlled way to a bosonic 
environment which has very similar properties to the radiation field within a photonic 
crystal. Here we study how, under certain conditions, and always within the two- 
level approximation, a controlled coupling with the environment can give rise to an 
improvement of the quantum search algorithm. To study the generality of this result, 
we analyze a real two-level system (as opposed to an effective one arising from the 
interplay of many qubits) which undergoes an adiabatic evolution, not necessarily 
corresponding to the search algorithm, and which interacts with the same photonic- 
crystal-like environment. We show that, even in this situation, the coupling with the 
environment makes the system, in certain parameter ranges, end up in a final state 
that corresponds, with a higher probability, to the target state (the ground state of 
the final Hamiltonian) . 

In both cases, the quantum open system is described through a master equation 
formalism. Hence, the effect of the environment in the system dynamics is encoded 
in a collection of terms that depend on the so called dissipation rates. In this paper 
we stress the importance of considering these rates as complex quantities. Indeed, 
the conclusions about the performance of the adiabatic algorithm highly depend on 
whether the imaginary parts are neglected or not, which suggests that they cannot be 
considered simply as an overall energy shift. We use units in which h = 1 . 

2. Adiabatic Grover's algorithm 

In Grover's problem, the goal is to find a particular item, m, in an unstructured 
database of size = 2". The best classical algorithm for this search takes time 0{N), 
while it was shown in |12| that there is a quantum algorithm solving the problem in 
time 0{\/N), known to be the optimal performance [TTl 118) . 

The quantum algorithm maps each element of the database onto one element of 
the computational basis for n qubits. The solution will then correspond to a particular 
state, |m). In the adiabatic algorithm |13j . the system is started in an easy to prepare 
initial state, containing an equal superposition of all basis states. 



where \z) are the elements of the computational basis. Then a time dependent 
Hamiltonian is applied which smoothly interpolates between the initial Hamiltonian 




(1) 



i?o = l-|*o>(*o|, 



(2) 
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having |5'o) as its ground state, and the final one, 

H,n = t-\m){m\, (3) 

whose ground state is the solution |m). If the adiabatic condition is satisfied, the final 
state will be j^*/) w |m). 

The Hamiltonian governing the evolution can be written as 

Hais) ^ (l - s)Ho + sHra, (4) 

in terms of the dimensionless parameter s € [0, 1], a function of time that controls 
how fast the Hamiltonian changes. Hq acts non-trivially only on the two dimensional 
subspace spanned by |m) and \m±) = Sx^m 1-^)' ^^'^ '^^^ ^^^^ diagonalized 

analytically. Indeed, in this subspace it can be expressed as 

HG{s) = h + His), (5) 

where we have defined H{s) = ^{n{s)cr^ - A{s)(t^) with A(s) = + (2s - 1) 

and n{s) = 2(s — 1) ^^^^ , ctx.z being the Pauli matrices in the basis |mj^)}. 

The time-dependent energy eigenvalues are Eq i{s) — T ^^^''^''2^'^^'''' ' ^^"^ 
corresponding eigenvectors 

|0(t)> = sin0|m) +cos0|m_L), 

|l(t)) = - cos^lrn) +sin6'|m_L), (6) 



where sin 9 — \/ a^d cos 9 = . . The subspace orthogonal 

V 2q(s) ^2a(s)(a(s)+A(s)) ^ 

to \m) and \mj_) is the eigenspace of Hq with energy E2 — I. The time dependent 
gap between the ground and first excited state is then given by a{s) = Ei — Eq = 

The adiabatic condition imposes a lower bound on the running time of the 

|0{s)) 

algorithm, T 3> max a(t)-^ ' '^here the maximum is taken over the whole range 

se[o,i]. 

If the interpolation is linear, i.e. sun{t) = t/T, and the total time is larger 
than Tiin — -^i the adiabatic condition is approximately satisfied, with an error 
in the overlap between the final state and the solution. Using a different function 
s(i), it is possible to ensure the condition by adjusting the velocity of change of the 
Hamiltonian to the instantaneous gap [13 . This leads to the optimal performance, 
when the interpolating function is 



1 L ^ tan(2e^^^ - arctan^/F^) ^ 
SoMt) - 2 (^1 + ^7^=^ j • (7) 

In that case, the total running time of the algorithm is Top* — tta/ZV /2e for ^ 1. 



3. Evolution equation of the system weakly coupled to an environment 

Let us consider the evolution of a system evolving adiabatically with Hamiltonian 
Hs{t), and interacting with an environment whose free Hamiltonian is Hb- Then, the 
total Hamiltonian has the form 



H — Hs{t) + Hb + -ffint' 



(8) 
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where H-^^^^ is the interaction Hamihonian, which describes a Unear couphng between 
system and environment operators. In particular, for an environment given as a set 
of harmonic oscillators, described by annihilation bx and creation 6^ operators, the 
interaction may have the form 

^int = EEffAMl+h.c.) (9) 

A i 

where g\ is the coupling constant of the system with the mode A, and ct; (aj) are spin 
ladder operators corresponding to the qubit i. 

In our case, we will consider that the environment is mainly coupled to the 
two-dimensional subspace spanned by |m) and \m±) [IHIII]. In this situation, the 
interaction Hamiltonian © can be written approximately [12] as ifjn^ = A ^ B, 
where A = az and B are operators that act on the system and the environment 
Hilbert spaces, respectively. Considering that the environment is composed of a 
set of harmonic oscillators, the environment coupling operator can be written as 
B = X^A ff^(^A + ^a)- The evolution equation of the system within the weak coupling 
approximation is given by (see the appendix for further details). 



dps 
dt 



- t[Hs{t), ps] - f drgit - r) (AA(r, -t)ps - A{t, -t)psA) 
Jo 

- [ dTg*it^ T)ipsAiT,~t)A~ ApsA{T,-t)) . (10) 
Jo 

with A{T,-t) = U{t)U^T)AU{T)U^t), U{t) rexp[-i Hs{T)dT] {T being the 
usual time-ordering operator), and 

9{t) = E^Afcoth (^^j cos{ujxt) -tsiniiuxt)]. (11) 

The well known Lindblad equation can be recaptured from the former equation by 
just considering a delta-correlated bath, so that g{t — t) = TS{t — r). In such a case, 
J^dTA{T,-t)g{t - t) = TA, and similarly drA^r, -t)g* {t - t) = T*A. Hence, 
the Lindblad equation is not valid for every system that is weakly coupled with its 
environment. It requires, in addition, that the coupling operators of the system (in 
our case there is a single one. A) evolve very slowly with the system Hamiltonian, 
within the time scale tq in which the correlation function g^r) decays. 



3.1. Bloch-Redfield equation 

The Bloch-Redfield equation describes the evolution of the density operator in the 
diagonal basis \n) of Hs, Pnm = (nT-lpsln), with pmn = {'m'\ps\n). In our case, the 
system Hamiltonian is time dependent, and the time dependent eigenbasis \n{t)) 
(corresponding to the set of eigenvalues En{t)) that diagonalizes instantaneously 
the system Hamiltonian Hs{t), should be considered. In this situation, the matrix 
elements of the system density operator are pmii — (w(i)|ps|"-(i)), and their evolution 
equation can be written as 

Pnm = {Tn(t)\ps\n{t)) + (m(t)\ps\n(t)) + (rn{t)\ps\h(t)) . (12) 

The first term on the right hand side of the equation is just the usual master equation 
([10]), expressed in the instantaneous basis of Hs{t). In order to project pO)) on the 
system basis, one should calculate terms such as A{T,—t)\n{t)), which requires the 
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calculation of quantities such as U{t)\n{t)) . Considering as a good approximation 
that the free system is undergoing an exact adiabatic evolution, so that the adiabatic 
theorem can be applied, we can write |20) 

U{t)\n{t)) =e-'/o^"(")''"|n(t)). (13) 
Hence the first term of the evolution equation (|12l) can be expressed as 

{n{t)\ps\p{t)) = - iEnpPnp + ^Tqn{t)AnqPqlAip + '^T*gAnlPlqAqp 

ql ql 

- ^ TqlAnlAlqPqp ~ ^ TpgAnlPlqAqp, (14) 
Iq ql 

where Anm = {n{t)\A\m{t)) , E^p ^ E„ - Ep, and 

Tqi{t)^ rdT9(t-r)e^/>"'(^'("')-^'("'» (15) 

"'0 

are the dissipation rates mentioned in the introduction. Notice that in our two level 
system, any energy difference is Eq{t) — Ei{t) — ±a{t), where q,l = 0,1. In order 
to further simplify the equations we may assume that the rate of variation of a{t) is 
much slower than the decaying of the correlation function g{t). In other words, when 
the condition 

1 da 1 

< — (16) 



a{t) dt Tc 

is satisfied, the integrand Eq{T') — Ei{t') — ±q;(t') appearing in (fTS)) will variate very 
slowly in the integration region, and we can write 

Tqi{t) = /* drgit - r)e'(^^(*)-^'(*))('-). (17) 
Jo 

Let us consider the evolution equations for our case, where the coupling operator 
A = az- Expressing the density operator ps in the spin basis {(Jx{t),(Jy{t),az{t)} as 
ps = ^(1 + PxO'xit) + PyCy{t) + PzCTzit)), thc master equation can be written as follows 



2csr^ - 4c^Re(roo)px + ait)py + 2 ( csr+ + i 
2cs( r+-2Re(roo) )pz 



Py = 2cs(^r+ - 2Im(roo)j + \^2s'Tj - a{t) )p, 
- 2 f 2c2Re(roo) + ^^T+Apy + 2csr7p, 



2-p+ 



Pz - 2s2r^ - 2(2csRe(roo) -9) + 2s'Rc{Tno))px + 2cs(r+ 

- 2Re(roo))p^, (18) 

where = Re(rio - Tqi), = Re(rio + Tqi), r7 = Im(rio - Tqi) and 
— Re(rio + Fqi). The dissipation rates are defined as 

Too = / dTg{T), 
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dTg{T)e 



— ia{t)T 



(19) 







and the variables c = cos 26* — — — , and s = sin 20 ~ — — . 
4. Adiabatic evolution in different environments 

We now analyze two different cases in which adiabatic evolution occurs in the presence 
of dissipation. 

In the first place we consider the effect of a thermal environment on the 
performance of the adiabatic search algorithm. A thermal environment is the one 
that we would expect to find naturally when a system is coupled to a bosonic bath, 
for instance phonons in a solid lattice, or the radiation field at finite temperature. 

In the second case we consider our system to be coupled in a controlled way 
to a bosonic environment that corresponds to a matter wave field. This example 
is much more specific than the former, but allows us to illustrate a case in which, 
provided that the two-level approximation remains valid, an artificially controlled 
coupling allows in some regimes an improvement on the performance of an adiabatic 
quantum computation. 

In all our simulations we consider the number of qubits n = 10. 

Thermal environment 

As noted above, the most important quantity that characterizes the influence of an 
environment on the system dynamics is the so-called correlation function. When the 
exact form of the coupling constants or the dispersion relation of the environment 
a;(k) is not known, a phenomenological model should be used to describe the 
interaction. In that situation, we can express the correlation function (|lip as 



which fulfills the property g{—T) — g*{T). In the last expression, the function J{uj) 
is the so-called spectral density of the bath. A very well known approach consists in 
assuming that J(w) behaves as uj^ [T5]. In that case, it can be written as 



This model of spectral function has been extensively studied in the context of the spin- 
boson model |15j . where three different regimes were described: a sub-ohmic regime in 
which < s < 1, a super-ohmic regime with s > 1, and an ohmic regime where s = 1. 
The exponential factor appearing in the last expression has been added to provide a 
smooth cut-off for the spectral density, which is modulated by the frequency ujc- This 
parameter controls the correlation time of the environment, approximately given by 
Tc ~ l/wci the larger lUc, the smoother the spectral function, and the shorter the 
time the environment takes to relax to equilibrium, giving rise to a more Markovian 
interaction. As seen in |15) . whether we have ohmic, sub-ohmic or a super-ohmic 
spectral function depends on the type of reservoir, and determines quite strongly the 
evolution behavior of the coupled system. In this work we will focus on one of the 
most significant cases, the ohmic dissipation, in which s = 1, and on large uJc, when 
the spectral function is a smooth function of the frequency. 




(20) 



J{uj) = riuj'ujl-'e-'^''^-' 



(21) 
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(a) Fig l.a. 



(b) Fig. l.b. 



Figure 1. Considering the coupling with a thermal environment, final success 
probability of the algorithm with adiabatic final adiabatic time T, as a function 
of T/Tmax (being Tmax = 0.8Tji„, with Tji„ = the time required for final 
success probability 0.55 in the closed case) for different coupling parameters rj, and 
for two cases: Fig l.a. considering complex dissipation rates Tij (with i,j = 0, 1), 
as they appear naturally in the derived Bloch-Rodfield equations, and Fig l.b. 
considering only the real part of the dissipation rates ^{Fij} (discontinuous 
lines). In both figures, the black solid line corresponds to the closed case rj = 0, 
while black, green and blue discontinuous lines corresponds to r; = 0.05,0.1,0.5 
respectively. The curve r] = 0.5 is not shown in Fig l.a., since for that value the 
imaginary part of the rates is too large and gives rise to unphysical results, what 
points out a failure of the weak coupling assumption. The frequency cut is chosen 
as uJc = 0.25. 



Let us study how the final success probabihty of the adiabatic search algorithm 
is affected by the presence of a thermal environment. To this order, we consider 
the evolution equations derived in the former section, with dissipation rates that 
depend on the correlation function (PD|) . As shown in figure l(a)[ the final success 



probability decreases in the presence of dissipation for any value of evolution time 
of the adiabatic algorithm T considered. A common approximation in the literature 
consists in neglecting the effect of the imaginary parts of the dissipation rates (fTO)) . 
with the argument that they amount only to an energy shift which in the context of 
quantum optics is known as the Lamb shift |21| . However, in the current scenario, 



figure 1 (b) shows that if only the real part of the dissipation rates is considered in the 
equations, the result for the same couplings is completely different. This shows that, 
at least in the present case, the imaginary parts of the rates cannot be considered as 
a simple energy shift. In fact, when eliminating the imaginary parts, one can see that 
the larger the coupling with the environment, the higher the final success probability 
at any adiabatic evolution time. 

As a consequence, we stress that in order to correctly describe the system, 
the imaginary parts of the rates should in general be considered, since they appear 
naturally in the derivation of the second order evolution equation of the system (in 
other words, they correspond to second order terms, as well as the real parts), and 
they produce relevant changes in the evolution. 

4-2. Structured environment 

The effect of dissipation on the performance of the adiabatic search will depend on the 
details of the environment and its interaction with the system. The particular features 
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of the adiabatic search make it specially sensitive to energy scales in the bath that are 
comparable to the characteristic gap in the system, AE [31 [TU] . Intuitively, one would 
expect that the algorithm is more resilient when this energy corresponds to a gap in 
the spectral function of the environment (i.e. a frequency region where this function 
is zero) , similar to the effect when the dominant energy scale of the bath is far away 
from AE', as studied for the thermal case in [1^. However, while this picture gives a 
qualitative image, a complete idea of the effect of the environment in the algorithm 
can only be achieved with a more complete study. With this aim, in the following 
we analyse in detail the effect of an structured environment in the performance of 
the search algorithm. Particularly, we will consider an environment with a gap in the 
density of states (which leads to a gap in the spectral density), and which has the same 
characteristics as the radiation field within a photonic crystal. In addition to that, 
it is specially interesting to consider the case of a controllable environment, where 
the capability to tune the parameters of the interaction may give us the possibility 
to improve the performance of the algorithm via the dissipative dynamics. Such an 
environment is indeed realizable in the context of optical lattices, as we discuss in the 
following. 

Let us therefore consider an ideal situation in which the quantum search algorithm 
is implemented in an atom lattice, which is formed by loading a gas of ultracold atoms 
in a standing wave field. In the last few years, atom lattices have been experimentally 
realized by several groups [22l [23l [24l [25l [26l [27] , and due to their high controllability, 
they have also been proposed as a candidate to realize a quantum computer [27j . 

In addition, it has been recently shown [28l [16] that atoms in an optical lattice 
can be coupled to an environment in a controlled way. In the later proposal, atoms 
are considered to have two relevant internal states a and b. Atoms in state a are 
actually trapped by the optical lattice, and have a frequency = + a;o/2, with 
Wo the trap frequency of the lattice. In addition, they are considered to be in the 
so-called strongly correlated regime, where either there is one or zero atoms at each 
site of the lattice. Atoms in state b are not trapped by the lattice potential, and 
have an energy huji, + where huJb and correspond to the internal and kinetic 
energies, respectively. If a Raman transition of total frequency uj^ and Rabi frequency 
fii is produced between the trapped and the untrapped state, the Hamiltonian that 
describes this process can be written in the interaction picture as 

^int = EEsfc (bla.e^^'^'-^^^-^^y'^^+h.c) . (22) 

j k 

In ([m, the sum in j runs over the N sites in the lattice, rj denotes the positions in 
the lattice, and Ak = /2m — A^,, with A^ = lul — {oJb ~ ^a) the laser detuning. The 
coupling constants are gk — riie"^"'' /'^{Sn^^'^Xj^/V)^^'^, where ~ (1/2towo)^/^ is 
the size of the wave function at each site. 

The previous Hamiltonian is very similar to the one describing the interaction 
of two level atoms with the radiation field. However, here the spin operators Ui 
(equivalently aj) are not describing transitions between atomic internal states, but 
rather they describe transitions from a Fock state (which corresponds to the 
presence of an atom at site i), to the Fock state |0)i, describing the absence of atoms 
at site i. On the other hand, the operators 6j^. (equivalently 6k) correspond to the 
creation (annihilation) operators of a bath of harmonic oscillators, which in this case 
is not the radiation field, but the matter-wave field that describes the untrapped 
atoms. 
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In our model, we assume that the qubits undergoing adiabatic evolution are 
encoded in the presence or absence |0)j of an atom at site i of the lattice. The 
coupling with the environment produces transitions between these basis states through 
the Hamiltonian (j22[) . However, contrary to the case of the thermal environment 
studied in the former section, here the coupling with an environment is not an 
uncontrolled event, but rather it is artificially produced through the two-photon 
Raman transition to the untrapped level. For that reason, several external parameters 
that describe the interaction can be easily controlled: most importantly JIl, which 
determines the coupling strength, , which determines the resonance condition, and 
Wo, which determines the trapping frequency of each lattice well, and as we will later 
see, will also characterize the width of the spectral function in the frequency space 

It is important to notice that, contrary to the thermal case, in this situation the 
correlation function can be fully determined by the coupling parameters gk and the 

.2 

dispersion relation of the environment, given by ujk = ^ + ^b- No phenomenological 
model is needed to characterize the spectral function. On the other hand, the relation 
uJk indeed resembles that of the radiation field in a three dimensional and infinite 
photonic crystal near the band-gap edge, which here corresponds to the frequency uJb 
|29[ 1301 131j . This environment will give rise to a very particular dissipation in our 
system. The correlation function can be written as [16] 

sW^^sie-^-'^f^i^, (23) 

k ^* 

where we have assumed that the sum in the wave vector k can be performed in the 
continuum limit. Here, the quantity vt = + i^ot. 

In figure 2(a) we observe, always in the two- level approximation and for various 
parameter regimes, a larger final probability of success at small values of T for the open 



system (discontinuous lines) than for the closed one (solid line). Similarly, figure 3(a) 



shows an improvement in the algorithm performance for certain values of A^. This 
results suggests that, at least within the two level approximation, a controlled coupling 
with a certain environment may improve the performance of the quantum search 
algorithm. Figures 2(b) and |3(b)] show the result for the same parameter regimes as 



2(a) and 3(a) respectively, but just considering in the equations the real part of 



the dissipation rates. Just as in the case of the thermal environment, we can see that 
results for real and complex rates differ considerably. 

Indeed, considering the full complex rates is particularly important for a 
structured environment like the one analysed here. When the system frequency is 
within the gap, the corresponding rate is, at long times, a purely imaginary quantity. 
Hence, by only considering the real part of the rates one would arrive to the wrong 
conclusion that the environment has no effect at all in the system. However, it 
has been known for a long time that even when the system resonant frequency is 
within a gap of the spectral density, the coupling with the environment has important 
consequences in its dynamics. Among other things, a so-called photon-atom bound 
state is formed [29] , in which the energy is coherently interchanged between system and 
environment. This particular state, and more generally the dynamics of the system 
with a resonant frequency within the environment gap, can only be properly described 
when considering the imaginary part of the rates. 

In addition we note that, while in the thermal case the improvement of the 
algorithm performance observed for real rates can no longer be observed when 
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(a) (b) 

Figure 2. For a system coupled with a structured environment, ttie figures 
siiow the final success probability of the algorithm as a function of T/Tmax (with 
Tmax = Tii„, and Tii„ = ^^ )- Solid black lines represent the solution for a closed 
system; red and blue dashed lines correspond to = 0.05 with = 0.2 and 0.27 
respectively; red and blue dotted lines correspond to ?7 = 0.1 with Aj^ = 0.2 
and 0.27 respectively. Fig l.a. represents the curves for complex rates, and for 
different couplings 7 and laser detunings A^. Fig 2. a. represents the curves for 
the same parameters but considering only the real part of the dissipation rates. 
In both cases the trap frequency is chosen as ujq = 0.25. 



(a) (b) 

Figure 3. For a system coupled with a structured environment, the final success 
probability of the algorithm is plotted as a function of the laser detuning A^. 
Fig. 3. a represents the solutions considering complex rates. In this figure, the 
solid black line represents the solution for the closed system, while black, green, 
red and blue discontinuous lines correspond to couplings rj = 0.01,0.05,0.1,0.4 
(from the bottom up) respectively. Fig. 3.b. shows the curves considering only 
the real part of the dissipation rates. In this figure, the solid black line represents 
the solution for the closed system, while black, green, blue discontinuous lines 
correspond to couplings -q = 0.01,0.05,0.2 (from the bottom up) respectively. 
The trap frequency is chosen as = 0.25, and Tmax = 0.8 Tun, with T;i„ = 



considering the full complex rates, the opposite is observed here. Indeed, as noted 
above only when complex rates are included in the description, an increase in the final 
success probability is observed for certain parameter regimes. This can be seen again 
by comparing figures 2(b) and |3(b)[ corresponding to real rates, with 2(a) and 3(a) 



From all this we again conclude how crucial it is to consider the imaginary terms 
of the dissipation rates in order to correctly describe the system dynamics. 
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4.3. Application to a genuine two-level system 

In the former section we have seen within the two-level approximation, how an 
adiabatic search algorithm can be improved by connecting the system to a photonic 
crystal like environment as the one described in [16]. We are aware that the 
applicability of this result to a future implementation of the adiabatic version of the 
Grover's algorithm depends on how good the two-level approximation is. The validity 
of the two-level approximation will be analyzed in more detail in Sect. (5] Notice, 
however, that this caveat would not apply to a genuine two-level system controUably 
coupled to a dissipative environment. Therefore in this section we will study if a 
real two level system undergoing a more general adiabatic evolution (not necessarily 
corresponding to the initial conditions of Grover's algorithm) would also experience 
an improvement in the final result under some condition. We recall that here an 
improvement means that at the end of the adiabatic evolution, the population of the 
ground state of the final Hamiltonian is larger than in the closed system case. 

To this order we consider the Hamiltonian (|22]) but with n — 1, corresponding to 
a single site in the lattice. This model can be realized by considering a lattice with a 
filling factor sufficiently low so that the sparse atoms do not interact with each other. 
Hence, in a real setup we would have several independent copies of atoms at single 
sites. Instead of the adiabatic search algorithm, which should be realized with N 
sites, let us consider an adiabatic process described by the Hamiltonian (|4]), but with 
|m) = |1) and I-^q) = ao|0) -I- 6o|l) with arbitrary coefficients oq and 60 = ^1 — Oq. 
We observe in Figure |3] that, for certain parameter regimes the final probability of 




,51 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

T/Tmax 



Figure 4. Final success probability (population of the ground state of 
the final Hamiltonian) for a genuine two level system undergoing adiabatic 
evolution, considering different adiabatic evolution times T, with Tmax = 10 (the 
approximate time in which the closed system reaches a success probability equal to 
one). The solid line corresponds to the closed system, while black, green and red 
discontinuous lines correspond to the results for rj = 0.01,0.05,0.2 respectively. 
Other parameters are chosen as luq = 0.5, A/, = 0.5, and ao = \/0.5. 

success for the open system is larger than for the closed system. Similar results have 
been observed for other choices of the coefficient oq. Hence, we conclude that the 
improvement of the adiabatic evolution for the system coupled with the environment 
is a robust effect in the two-level system, not linked to the particular choice of 
initial state in the search algorithm. Although in this section we have studied in 
particular a system subject to a photonic crystal-like environment, we believe that 
similar results could be encountered in other different scenarios, provided that the 
interaction between the system and the environment can be externally controlled. A 
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simple example would be an atom in which the transition between the two relevant 
internal levels is produced by a two-photon process, mediated by a single laser followed 
by a spontaneous emission. This corresponds in quantum optics to a A scheme, and 
would give rise to an effective Hamiltonian of the form (1^^ . Despite having a different 
dispersion relation, this Hamiltonian would describe a coupling with the environment 
that can be tuned through the laser Rabi frequency. 

5. Spectral density and the validity of the two level approximation 

Let us now qualitatively study the validity of the two-level approximation. This 
approximation is suitable when the environment does not produce significant 
transitions to other levels of the energy spectrum different than |0) and |1), the ones 
involved in the quantum search algorithm. In general, the transition probability per 
unit time between two system levels \i) and \j) is approximately given by the Fermi 
Golden Rule as Pij « J{uJij){i\H^^^\j). From this formula, we can already see that 
the spectral density, is indeed a very important quantity that determines, at 

each frequency, what is the density of environmental states available, and how strong 
the coupling is. Indeed, one may expect that the two-level approximation applies 
better to parameter regimes such that the spectral density is small for frequencies 
corresponding to transitions to levels different from |0) and |1). 

Notice that this can only be understood as a first qualitative approach to the 
question of the validity of this approximation. Firstly, because an equally important 
factor to determine Pij is the magnitude of the transition amplitudes {i\H[iif^\j) ^ and 
secondly, because Pij only accounts for the real part of the transition rates. 

Let us illustrate this qualitative approach by considering the energy vs. time plot 
in figure [5] for the thermal bath. All energy differences Ei — Ej in the closed system 
are represented in the plot as solid lines. The background shows how the spectral 
density varies for different energies, with darker color representing a larger value of 
J. Intuitively we expect that a particular transition i ^ j becomes more important 
when the corresponding energy gap coincides with large values of J{i-o)- 




0.2 0.4 0.6 0.8 



Figure 5. The energy gaps in tlie closed system as a function of the adiabatic 
parameter s: E\ — Eq (blue), E2 — E\ (red) and E2 — Eq (green). For each value 
of energy, the background color indicates the magnitude of the spectral density 
for a thermal bath, where we have chosen the same parameters as in section l|4.1| l. 
s = 1 and LJc = 0.25. 

In order to compute the spectral function for the photonic crystal-like 
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environment, let us consider the calculus of the correlation function of the system 
as defined in (j^ . 

Here, both and cjk only depend on the modulus of the wave vector. Using this, 
the angular part of the integral has been performed (giving rise to a factor 47r) and 
only the integral of the modulus is left. We can consider that the coupling parameter 
does only depend on the modulus of the wave vector, because it has been assumed 
that the total wave vector of the Raman lasers sa 0. According to the definition 
(|20p. the spectral function is the kernel of the integral ([M)) translated to the frequency 
space. Hence, it is necessary to use the dispersion relation oj{k) — ^ + Wf, in order 
to perform a change of variable in the integral (1^^ , such that it becomes an integral 
in w, 



g{t) = T] doj^/2{oj - Wb)e-'("'-"'^)*e"^^^, (25) 



g^l/2Q2 

with ry = — 372^- Considering now an energy shift of the form uj — uj^ujl, the former 
integral is just 



g{t) = Tj dioy^2{uj - AL)e-*"*e~^^^ . (26) 

Since in this case the temperature of the reservoir is zero, it is straightforward to see 

that the spectral density has the form J{uj) = riy/2{uj — Ai)e '^0 . 

Like in the thermal case, figure [6] represents all energy differences Ei — Ej for the 
closed system in solid lines, compared to the spectral density in the background. From 
a naive interpretation of the picture, one would consider that, since the minimum gap 
El — Eq corresponds to an energy within the environment gap, where the spectral 
density is zero, the effects of the environment are somehow minimized. However, 
as it is shown in figure 2(a)[ for the same A/, = 0.28, the performance of the 



algorithm is even improved with respect to that of the closed system. Hence, the 
environment gap does not protect the algorithm from the interaction, but rather it 
produces some positive effects on its performance. As already mentioned, another 
important factor that determines the transition probability to other levels \i) different 
to the ones considered in the two level approximation (|0) and |1)) is the amplitude 
of the transition elements {i\H^^^\0), (z|i7jj^^|l), and in general {i\H^^-^\j), with H:^j^^ 
given by ((^ . However, this analysis is out of the scope of the present paper, and will 
be made elsewere. 



6. Conclusions 



We have studied the performance of the adiabatic quantum search in a dissipative 
environment. Our derivation of the Bloch-Redfield equation allows for a more complete 
account of the effects of the complex dissipation rates than previous studies. In 
particular, we have shown that neglecting the imaginary parts of such rates give rise to 
a completely different dynamic than when considering the full complex rates. Hence, 
conclusions about whether there is an improvement of the performance of the adiabatic 
algorithm within an open system, cannot be extracted by a partial analysis that only 
accounts for the real part of the dissipation rates. 
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Figure 6. The energy gaps in the closed system as a function of the adiabatic 
parameter s: Ei — Eg (blue), E2 — E1 (red) and E2 — E0 (green). The background 
indicates the magnitude of the spectral density for the photonic crystal-like 
environment, with Aj^ = 0.28 and ljq = 0.25 as typical parameters used in 
section 1 14.21 1. We expect that the dominant transitions induced by the couplig 
to the environment are <-> 1. The transition 1 <-> 2 may also contribute, but 
it will be suppressed during the first part of the algorithm because the Ei level 
is initially not populated. It will again be important after s 0.65. With these 
parameters, moreover, the minimum gap, around s = 0.5, is decoupled from the 
frequencies of the bath, what can be expected to bring some protection to the 
most sensitive part of the algorithm. 

Particularly, we have analyzed the effect of the dissipation on the adiabatic 
version of Grover's algorithm for two different settings, namely a thermal bath and a 
controllable environment with tunable parameters. While for a thermal environment 
no improvement is observed in the performance of the algorithm, coupling the system 
to a structured environment gives rise to a final success probability that for certain 
parameter regimes is higher than in the closed system case. Hence, we have found an 
example in which the performance of the search can be improved with respect to the 
closed system, by tuning the bath parameters appropriately. 

Our study, valid for the quantum search in the framework of the two-level 
approximation, can be extended to the adiabatic evolution of a genuine two- 
dimensional system. 
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of Excellence MAP, the SFB 631, and the DFG Forschergruppe 635, and by the EU 
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7. Appendix A 

In this appendix we proceed to derive equation (jlOp in the paper. The total density 
matrix can be written in the form 

Ptot (*) = PS (t) (g> Pb {t) + Pcorrel' (2 7) 

where ps = TrBipiQ^it)) and pB = Trs{piQ^(t)) are the reduced density matrices of 
the system and the bath, respectively, and P^orrel describes the correlations existing 
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between system and environment. We assume that in general we are dealing with a 
thermal bath, so that pb — e^*^^^, with /3 — 1/{KbT), Kb the Boltzmann constant 
and T the temperature of the reservoir. We consider the master equation for the 
reduced density operator in the interaction picture with respect to both the system 
and the environment, psi — U\t)pslA{t), where lA{t) is given by the time-ordered 

exponential = Te~'/o Considering that = iHs{t)W(t), the 

evolution equation up to second order in g = 1 1 ^fjjit 1 1 / 1 1 -^-S + Hb\\ is given by [35] 

psi{t) = - I dTTrB[Hi{t), [Hj{T),psi{t) ® pb]]- (28) 
Jo 

Here the terms ^ TrB{pBHf_^^^) have been neglected, an assumption that is valid for 
most types of environments and couplings. Particularly, it is valid for our present case 
of a system lineraly coupled to an environment of harmonic oscillators. In addition, 
it has also been assumed that p{t) — ps{t) ® Pb, and since tc <C l/F, we have 
also assumed that psi{t) « psiif) in the integral. On the other hand, the quantity 
Hi{t) = W{t)Hl^^U(t) with Hf^^ = e*-f^«*iJinte"*^''*- Going back to the original 
picture, the evolution equation of ps is then given by 

^= ^^[Hsit),ps] 

dTU{t)TrB[Hj{t), [Hi{T),psi{t)®pB]]U\t). (29) 



The latter equation can also be expressed as 

^ = ~i[Hs{t),ps] - ^ drTrB [hI^^, [Hiir, ~t), psit) (g> pb]] ,(30) 

where Hj{t, —t) = U{t)U^ {t)HT^^{t)U^ {t) . Expanding the commutator on the right 
hand side of the last equation, we find the following terms, 

Ttb [Hi{T,-t),ps{t)pB]]) = TrB{{Hl^^Hi{T,-t),ps®PB} 

- Hl^^ps®pBHi{T,-t)-Hi{T,-t)ps®PBHl^^). (31) 

Let us now consider our case, in which the interaction Hamiltonian can be written 
as -ffjj^^ = A® where A — and B ~ X)a •'?^(^a ^^)- '^^'^ make further 
simplifications in equation pip 

TrB{[Hl^^AHi{T,~t),ps{t)®pB]]) 

- g{t - r) {AA{t, -t)ps - A(t, -t)psA) 

+ g{T-t){psA{T,-t)A-ApsA{T,-t)), (32) 

where A{t - t) ^ U{t)U^ {T)AU{T)UHt), and g{t - r) = TrB[B*B^pB], g{T - t) = 
TrB[B'^ B* Pb], having the form 

git) = ^ 5a [coth (^^^ cosiujxt) - isin(wAi)]. (33) 

Here, we have taken into account that, for a thermal bath, 
TrB{pBb{by) = Sx^x'N{ujx), 

TrsipBbxbl,) = 1 + h,x'N{ux), (34) 
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with N{ll!x) — l/(e"*'^ — 1) the number of excitations with frequency ujx- Note that 
according to equation ([Tl]) . g{~t) — g*{t). Considering that, and inserting in 
we find 

^ = - i[Hs{t), ps] - drgit - t) {AA{t, -t)ps - A{t, -t)psA) 

- f dTg*{t~T){psA{T,~t)A~ApsAiT,~t)). (35) 
Jo 
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